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We consider both the effect of particle inertia on stochastic Stokes' drift, and also 
a related process which could be considered as a crude model of stochastic Stokes' 
drift driven by an eddy diffusivity. In the latter, the stochastic forcing is a stable 
Ornstein-Uhlenbeck process rather than Brownian motion. We show that the eddy 
Stokes' drift velocity has a peak at a non-zero value of the correlation time-scale for 
particles that have the same (limiting) diffusivity. For both of the models considered, 
this study shows that not only can stochastic Stokes' drift be used to sort particles 
with different diffusivities, but also it can be used to sort particles of the same 
diffusivities but with different particle masses or correlation time-scales. This effect 
may be important in particle sorting applications. 
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1. Introduction 

This study is a continuation of the work of Jansons & Lythe (1998) on stochastic 
Stokes' drift, which is the modification of the classical Stokes' drift of a particle 
in a travelling wave due to the effect of Brownian motion (or some other random 
forcing with very short timescale correlations). In Jansons & Lythe (1998), the 
particle inertia was neglected. In this study, we first consider, in §2, the effect 
of particle inertia on the stochastic Stokes' drift velocity. For simplicity, we still 
neglect the inertia in the fluid, so we can use the Stokes flow approximation to 
determine the drag on the particle, and so that this drag is not history dependent. 
This is a consistent approximation for a particle in a gas, where it is reasonable to 
assume that the particle density is much greater than that of the ambient fluid. For 
the motion of a particle in a liquid of similar density to that of the particle, the 
situation is much more complex, and the work here is not applicable, though may 
agree qualitatively in some cases. 

Many authors have continued the investigation of various aspects of stochastic 
Stokes' drift. The first exact solution to Brownian motion forced Stokes' drift is 
due to Van Den Broeck (1999), who found an exact solution for an arbitrary wave 
in an integral form and simple exact solutions for a square wave and an impulse 
wave. Applications of stochastic Stokes' drift to ocean circulation have recently been 
considered by Restrepo (2006), and the techniques considered in this study appear 
to be relevant. Li Yu-Xiao et. al. (2001) studied the effect of asymmetric potentials 
on the stochastic Stokes' drift. Bena et. al. (2005) considered, both analytically and 
numerically, Stokes drift forced by a dichotomous Markov process. 
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Also the theory of stochastic Stokes' drift applies to Brownian motors and tilting 
ratchets (see Reimann (2002)). 

Jansons & Lythe (1998) showed that particles of different diffusivitics could be 
separated by the superposition of several waves with different directions, wavenum- 
bers and frequencies. In such a system, it is possible to arrange for particles to drift 
in different directions depending on their diffusivities. We show here that it is also 
possible to separate particles of the same diffusivity but with different masses. It 
is not yet clear if this effect will have practical applications, but the early signs are 
promising. 

In §3, we consider a related process of independent interest, and which could 
be considered a crude model of a particle moving in a travelling wave and small 
eddies, i.e. the random motion has short time-scale correlations. This is mathe- 
matically stochastic Stokes' drift forced not by Brownian motion but rather by a 
stable Ornstein-Uhlenbeck process. We show, for example, that it is possible for the 
stochastic Stokes' drift velocity to be increased by a careful choice of the correlation 
time-scale, even for particles with the same diffusivity in the long-time limit. Such 
an effect might be important in oceanographic applications of particle dispersion 
for example. 

In §4, we compare the leading-order asymptotic approximations of the Stokes' 
drift for the two types of models considered here with Monte Carlo simulations, 
which agree well even for reasonably large values of the 'small' asymptotic param- 
eter. 

Finally, in §5, we consider applications of these results to particle separation 
schemes using multiple waves in dimensions greater than one. 



As shown by Jansons & Lythe (1998), stochastic Stokes' drift in higher dimensions 
with many waves is, to leading order, a superposition of the contributions from 
each wave separately. This is because, in the leading-order calculation, cross-terms 
from different waves average to zero if they have either a different spatial or different 
temporal frequency. This is also true for the generalizations considered in this study. 
So, for simplicity, we consider the essentially one-dimensional system with a single 
wave, and comment on the higher-dimensional versions in §5. 

The standard form of the Langcvin equation for the velocity U of an isotropic 
Brownian particle in a gas at rest is given by 



where m and b are respectively the particle mass and particle mobility, both as- 
sumed constant, and i] is a time-dependent noise term. 

Note that assuming that & is a constant follows from the assumption that we 
can neglect inertia effects in the motion of the fluid around the particle. This is a 
rational approximation in the case of a particle in a fluid, provided that the particle 
density is much greater than that of the fluid, which is normally so for a gas, but 
is rarely true for a liquid. In the case of a liquid, the particle drag term depends on 
the history of the particle's motion (see Hinch (1975)), and we need to include both 
fluid and particle inertia terms in the description of the motion, which we shall not 
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do here. However, the results for a particle in a liquid are likely to be qualitatively 
similar to those of this study. 

Turning (2.1) into an ltd stochastic differential equation we find 

dU t = -XU t dt + \adB t , (2.2) 

where B is a standard Brownian motion (i.e. a Wiener process), A = (fern) -1 , 
a = {2bKT) 1 / 2 , K is Boltzmann's constant, and T is the (absolute) temperature. 
The long-time particle diffusivity is \u 2 . 

We now introduce an additional velocity term ef due to some wave motion 
imposed on the fluid, giving 

dU t = -X(U t - sf(X t , t))dt + XadB t , (2.3) 

where X t is the particle position and e is a dimensionless constant, which will be 
used as a small parameter in the asymptotic analysis below. The stochastic Stokes' 
drift velocity is defined as 

V = lim r 1 {X t - X ) . (2.4) 

t—>oo 

Using the Ergodic theorem we are able to replace this limit with an expectation of 
Ut, provided the starting phase of the wave is uniformly distributed. 

For a wave with characteristic wavcnumber k and angular frequency u, the 
natural time-scales in this system are 

A" 1 , lo-\ a- 2 k- 2 , (2.5) 

which are respectively the relaxation time-scale of the particle velocity, the period of 
the wave and the time for the particle to diffuse over a wavelength. These time-scales 
give two independent non-dimensional constants, in addition to efok/uj, where fo is 
a typical value of /. We assume that all of the time-scales in (2.5) arc comparable, 
and denote the common time-scale as T. 

We can expand both U and X as a formal asymptotic series in e, in the limit 
e — > 0, namely 

U t = U^ +eC/ t (1) +e 2 C/ t (2) + ••• (2.6) 

and 

X t - X[ 0) + eX^ + e 1 xf ) + ■■■ (2.7) 

These local approximations are valid for an o(e^ 1 T) time range, which is sufficient 
for the determination of stochastic Stokes' drift as the exponential tails beyond this 
time range do not contribute to the leading-order asymptotic results. At times of 
order e~ x T some contributions switch order, but we do not need to address this 
point as we are not going to determine the first correction to the stochastic Stokes' 
drift in this study. 

We now consider the formal e expansion an order at a time. 

(a) Orders 

A _1 dI7 t (0) = -U^dt + adBt, (2.8) 
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which is a stable Ornstein-Uhlenbeck process and can be solved exactly in the form 

U[ 0) = e- xt W (±A<7 2 e 2A <) , (2.9) 

where W is a standard Brownian motion, with 1^(0) = 0, which means that we 
have started with its stationary law. Note that has a zero mean. 
Thus the autocorrelation of is given by 

Ep^U^} = \\o 2 exp(-A|s - t\). (2.10) 
(6) Orders 1 

= -UP + f(X?\t). (2.11) 

C/ t (1) = A / cxp(-Ar)/(X t ( ° ) T , t - r)dr. (2.12) 
Jo 



Thus 



(c) Order e 2 

It is at this order we find the first contribution to the stochastic Stokes' drift. 

\- 1 ^ r = -UW + f { xP,t)X?\ (2.13) 
where the dash denotes differentiation with respect to the first argument. Thus 

r(2) 



This gives 



and so 



(2) 



/ cxpi-X^fiX^t-rjX^Ur. (2.14) 
Jo 



= ^ y cxp(-Xr)f(xP T ,t - r)uP a dadT 7 (2.15) 



A 2 



(2) ,*oc poo poo 



/ cxp(-A(r + P))f\x\%,t - T)f(X<® t -a- (3)d(3dadr. 
Jo 

(2.16) 

This is about as far as it is useful to go in general, so we now focus on the particular 

case 

f(x,t) = ucos(kx — wt + <j>), (2-17) 

where u, k, u) and <j> are constants. The phase <j) of the wave will not alter the 
long-time limits of interest, so it is not important. However, for convenience, and 
to avoid starting transients, we choose <j> to have a uniform distribution on [0, 2ir). 
Note that in this special case, £[£/ t (0) ] = £[[/ t (1) ] = 0. 

So the leading-order Stokes' drift is from U^ 2 \ From (2.16), we find 



(2) /.oo roc roc 

= — / / cxp(— A(r + f3)) sm((j> s ) cos((j) c )df3dadT, (2. 
Jo Jt Jo 



X 2 u 2 k 



18) 
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where 4> s — kxj:°^ T —Lu(t — T) + (/) and (j> c — kX^} oi _p — uj{t — a — (3) + <f>. On averaging 
we find 



(2)i poo poo poo 



\\ 2 v?k 



pOO pOO pOO 

= -/ / / exp(-A(r + /3))E[sm(<j) s - <j) c )]d(3dadT. (2.19) 
Jo Jt Jo 



This is exact, rather than a long-time limit, since we chose (j> to have a uniform 
distribution, so the time t = is not special. 
Note that 

4> s - <$> c = k{x[% - x[%_ p ) +oj(r-a- (3), (2.20) 

which is Gaussian, and so E[sm(<p s — <f> c )] can be evaluated exactly. From (2.9), we 
see that 

C(t -s) = E[(X? X°f] = o*(\t - s\ + (e " A '7'" 1) ) • (2-21) 

Thus 

E[sm((f> s - 4> c )\ = cxp(-i/c 2 C(r - a - 0)) sin(w(r -a- (3)). (2.22) 

Putting everything together, we find from (2.19) the stochastic Stokes drift velocity, 
V, to leading-order in e: 



V 



CO I'OO pOO 



\e 2 \ 2 u 2 k J J T J ^ ^' 2 



cxp(-A(r+/3)-i/c 2 C(a+/3-T)) sm(Lu(a+/3-T))dpdadT. 

(2.23) 



We may simplify this a little, to give 

/>OC />0O 

V=\e 1 \u 2 k\ / exp(-A/3- ±k 2 C(a + f3))sm(uj(a + f3))df3da. (2.24) 
Jo Jo 



We should expect (2.24) to reduce to the result of Jansons & Lythe (1998) in 
the limit A — ► oo, i.e. in the zero mass limit. Taking this limit reduces (2.24) to 

V=\e 1 v 2 k\ exp(-±k 2 C(a))sm(uja)da, (2.25) 
Jo 

and 

C(a) = cr 2 \a\. (2.26) 

Equations (2.25) and (2.26) agree exactly with the penultimate line in the calcula- 
tion of Jansons & Lythe (1998), but in a slightly different notation. 
To leading order 

hm^S=a 2 . (2.27) 

i^oo t 

We could easily determine the first correction to the variance of position, but we 
do not present this result here as it was felt that it was a little too messy to be of 
interest. 
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3. Eddy Stokes' drift 



In this section, we consider a closely-related process to the stochastic Stokes' drift 
with inertia process of §2. This process is stochastic Stokes' drift forced by a stable 
Ornstein-Uhlenbeck process rather than Brownian motion. The governing equation 
(using notation as close to §2 as possible) is 



dX t 
dt 



= U t + sf(X t ,t), 



(3.1) 



where Ut satisfies (2.2), but, of course, does not have the same physical interpre- 
tation. This process could be considered as a crude (or qualitative) description of 
stochastic Stokes' drift with the randomness coming from eddy diffusivity. 

In the case f(x,t) = ucos(kx — u)t + <f>), the stochastic Stokes' drift, to leading 
order, is given by 



V 



12 2 



k / sin(wt)exp (-\k 2 C{t)) dt, 
Jo 



(3.2) 



which follows from an argument that is a simpler version of that in §2, and is the 
same as that of Jansons & Lythe (1998) but with <r 2 t replaced by C(t), so will not 
be given here. In this case, it is also trival to see that the A — > oo limit of this result 
is exactly the same as that for Brownian forcing in Jansons & Lythe (1998). 

We now determine the leading-order correction to the long-time limit of the 
variance of position. This will appear at order e 2 . We take = for all n. Note 
that, to the required order, t~ 1 E[X t ] — > 0, as t — > oo, so we need only to determine 



E[X, 



E 



(0K2 



+ 2eE 



which we consider term by term. 



r x E 



as t 



as t 



oo; 



oo; 



t~ l E 



(0)n2 



2 (e [(Af y 

t^Cit) -^a 2 , 
0, 



2E 



x (0) x (2) 



(3.3) 
(3.4) 
(3.5) 



r x E 



(^ (1) ) 2 



as t — * oo; 



as t — > oo. 
Thus 



= t-V 



COS 

(kX^ - aja^j cos (kxf ] - Lof3"j 

/ cos(t^) cxp (-\k 2 C{t)) dt, 
Jo 



dadf3 



r x E 



x t° )x t 2) 



0, 



(3.6) 
(3.7) 



hm — - = a + e 

t^oo t 



/•OO 

u 2 \ cos(ut) exp (-\k 2 C(t)) dt + -- - . (3.8) 
Jo 



Article submitted to Royal Society 



Stochastic Stokes' drift with inertia 



7 



0.018 




Figure 1. Stochastic Stokes' drift velocity against A for e = 0.2, and all other parameters 
set equal to 1. The curve is the leading-order asymptotic result and the points are Monte 
Carlo results. 



We can also consider the limit A — > oo of this result, giving the corresponding 
result for Brownian forcing, though this was not determined by Jansons & Lythe 
(1998) in their study of this process. In this limit, we find 

,. Var[X t ] 2 2 2u 2 k 2 a 2 

lim — = ° +e 2 —r— A „+■■■. 3.9) 



4. Numerical results 

To determine how small e needs to be for the asymptotic results to be a good 
approximation to the exact results, we perform some Monte Carlo simulations. In 
the simulations shown here, we use the stochastic Euler method with dt — 0.001 
and all parameters other than e set equal to 1. These simulations are run to t = 10 8 
to control the variance sufficiently for comparison of the time-averaged stochastic 
Stokes' drift. 

We present numerical confirmation of the results of both §2 and §3. 

In the case of stochastic Stokes' drift with inertia, and for e = 0.2, the agreement 
is excellent and is shown in Figure 1 , though the (predicted) variance leads to a little 
scatter around the asymptotic result. This could have been reduced by running the 
Monte Carlo simulations up to a a larger time t, but as it was each point took 5 
hours on a 3GHz Intel CPU. 

When e = 0.5 for this model, the scatter of the Monte Carlo results is reduced 
(see Figure 2) and a small departure from the asymptotic approximation can be 
seen, but it is still small enough for the asymptotic result to be useful in applications. 



Article submitted to Royal Society 



8 



K. M. Jansons 



0.12 




Figure 2. Stochastic Stokes' drift velocity against A for e = 0.5, and all other parameters 
set equal to 1. The curve is the leading-order asymptotic result and the points are Monte 
Carlo results. 

The agreement between Monte Carlo results and asymptotic results is much the 
same for the case of eddy Stokes' drift (§3). However, for this model there is a peak 
in the Stokes' drift, which was not seen in the model of §2, and is shown in Figures 
3 and 4 for e = 0.2 and e = 0.5 respectively. Thus the eddy Stokes' drift at finite 
A can be larger than that of the A — > oo limit (i.e. larger than for the Brownian 
Stokes' drift result). This does not appear to occur in the model of §2. 

Finally, we plot the first correction to the large-time limit of the particle variance 
for eddy Stokes' drift in Figure 5. 

5. Multiple waves and higher dimensional results 

In the work of Jansons & Lythe (1998), they considered stochastic Stokes' drift in 
several dimensions. In this case, it was possible to sort particles by a combination 
of waves in different directions and different spatial or temporal frequencies. At 
leading order in the Stokes' drift, such combinations of waves act independently 
on the particle as the cross-terms average to zero in the long-time limit. Thus the 
modifications to the classical Stokes' drift due to diffusion is different for each of 
the wave components and so the linear combinations of these drifts will in general 
result in particles of different diffusivities diffusing in different directions. This leads 
to an interesting continuous particle sorting method. 

In the current study too, the higher dimensional versions work in exactly the 
same way, as for the same reason all cross-terms due to different waves average to 
zero, at leading order, provided that they have either different spatial or temporal 
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Figure 3. Eddy Stokes' drift velocity against A for e = 0.2, and all other parameters set 
equal to 1. The curve is the leading-order asymptotic result and the points are Monte 
Carlo results. 
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Figure 4. Eddy Stokes' drift velocity against A for e — 0.5, and all other parameters set 
equal to 1. The curve is the leading-order asymptotic result and the points are Monte 
Carlo results. 
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Figure 5. The first correction to the variance coefficient for the particie position for eddy 

Stokes' drift against A. 



frequencies. Thus again the resulting leading-order stochastic Stokes' drift is just a 
linear combination of the stochastic Stokes' drifts of each component. 

This implies that, in the case of stochastic Stokes' drift with inertia (see §2), it 
is possible to sort particles that have the same diffusivity but different masses, by 
arranging them to have different directions for the resulting Stokes' drift. Whether 
this has practical applications remains to be seen, but it certainly looks hopeful. 

Equally, in the case of eddy Stokes' drift (see §3), particles with the same dif- 
fusivity in the long-time limit but with different correlation times can be sorted by 
arranging the resulting Stokes' drift to be in different directions. 



6. Conclusions 

In this study, we have extended the work of Jansons & Lythe (1998) in two di- 
rections. Firstly in §2, we consider stochastic Stokes' drift including the effect of 
particle inertia, and secondly, we consider stochastic Stokes' drift forced by a stable 
Ornstein-Uhlenbeck process rather than Brownian motion. Both of these results are 
to leading-order in the strength e of the wave motion, and the stochastic Stokes' 
drift appears at order e 2 . These results extend trivially to higher dimensions and 
many waves, as the cross-terms for waves of different spatial or temporal frequen- 
cies average to zero. Thus the 'fanout' of particles of different diffusivitics observed 
by Jansons & Lythe (1998) extends to particles of the same (long-time limit) of 
diffusivity but with different masses in the case of the model in §2 and different 
correlation times in the case of the model of §3. In the case of eddy Stokes' drift (see 
§3), we also compute the long-time limit of the variance of particle positions. The 
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short correlation time limit of this result applies directly to the Brownian motion 
forcing case studied by Jansons & Lythe (1998). 

The agreement of these results with Monte Carlo simulations is impressive (see 
§4), and shows that, for practical applications, the leading-order asymptotic ap- 
proximations arc sufficient, even when the 'small' asymptotic parameter e = 0.5. 
Note that this is partly due to the asymptotic expansion being naturally in e 2 rather 
than e, as replacing e with — e is equivalent to a phase shift, and does not effect the 
stochastic Stokes' drift velocity. 

It appears that Van Den Brocck's (1999) exact method could be extended to 
consider the system of this study, but it is less clear if there are simple exact 
expressions even for a square wave. 

One extension of these results, which would be of interest, is to consider the 
behaviour of flexible particles that are large enough to feel different parts of the 
wave. The simplest such example would be a flexible dumbbell, consisting of two 
point particles connected by a linear spring. 
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